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0^ ■ Abstract 



We demonstrate that coulomb interaction between a static charge source 
and and K~ from (j) meson decay could introduce systematic shift in the (j) 
invariant mass reconstructed from the K^K^ pair. 



a^ : 1 Introduction 



Hot and dense matter is created in heavy ion colhsions at the BNL AGS and the 
CERN SPS. It is predicted that characteristics of hadrons, such as their masses and 
widths, can be modified in hot and/or dense matter |l|. Such modifications are of 
^ , importance because they signal chiral symmetry restoration predicted by QCD 0. 

2 ■ The (j) meson is an ideal candidate to look for such modifications for the following 

reasons. First, mesons are short-lived, therefore relatively large fraction of (/)'s decay 
^ ■ inside the charge fireball created in heavy ion collision; meanwhile, the mass peak 

H : 

. 5t , is still narrow enough to be sensitive to small shift in the mass. Second, the decay 

momentum of — K^K^ is very small, therefore momenta of the secondary kaons 
do not contribute significantly to the invariant mass, resulting in little sensibility of 
any experimental uncertainty in the kaon momentum measurement on the mass. 

The meson was measured by the BNL AGS experiment 859 in Si+Au collisions 
at 14.6 A-GeV/c An analysis of the data showed that the invariant mass 
from K^K^ pairs detected in the E859 spectrometer decreased systematically with 
increasing centrality, reaching a value about 2 MeV smaller than the data-book 
mass in the most central events 0. The decrease was even more pronounced for 0's 
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at low p±. No experimental systematics have been identified so far to be responsible 
for the decrease in the invariant mass. 

In this paper, we will demonstrate that electro-magnetic coulomb interaction be- 
tween the positive charge fireball and the and K~ from a (p decay could give 
rise of systematic shift in the (p invariant mass, which is on the same order of that 
observed in the experiment. We will first give the formalism in Sect. 0, and then 
describe in Sect. ^ the Monte Carlo procedure we used implementing the experimen- 
tal acceptance. We will present the results of our calculation in Sect. ^, followed by 
discussion of limitations of the analysis. In Sect. ^ we will give our prediction of 
(f) mass shift due to coulomb effect in central Au+Au collision, and finally we will 
conclude in Sect. ^. 

2 Relativistic Analysis 

Suppose a static, spherical and uniform charge fireball is created in central heavy 
ion collision. (We will present all the formula in this section in the rest frame of the 
fireball.) The coulomb potential produced by the fireball at a distance r from origin 
of the fireball is given by 



for positive unit charge particles, where Z and R are the charge and the radius of the 
fireball, and a = ^ is the electro-magnetic coupling constant. 

The coulomb force will put a ± impulse onto the momentum. Due to the 
momentum impulses, there will be a shift in the invariant mass, the square of 
which is given by 



where p+ and E^, p_ are energy and momentum- vector of the and K , 
respectively. Thus, we can obtain the mass shift. 
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where di and Of are the opening angles at the decay point and at oo where the kaons 
are measured. Expanding Eq. ^ to the first order of V{r) and assuming that the 
opening angle does not change and the kaon momenta are large, we obtain 
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2p+p_ \p- p+ J '2^p+P- 
where (5± = p±/ E± and mx is the rest mass of kaon. As seen from Eq. ^, the mass 
shift vanishes on the first order when the kaon momenta are equal, and is negative 
for p_ > p+ and positive for p+ > However, the second order effect is always 
positive. Therefore, averaging over all possibilities of the (f) decay kinematics, the 
invariant mass will have slightly positive shift. 

Consider a central Si+Au collision. A fireball with positive net charge consisting 
of almost all protons from the Si projectile and a considerable fraction of protons from 
the Au target is formed. Using the simple spectator-participant geometry model of 
heavy ion collision, we can easily calculate the number of target participant nucleons, 
hence the amount of charge in the fireball, Z = 43. Assuming that the nuclear density 
of the fireball at freeze-out is as same as the normal nuclear density, we obtain the 
radius of the fireball, R = 5.2 fm. Suppose the decays at origin and neglect change 
in the opening angle, then we can calculate the mass shift according to Eq. ^ as a 
function of p+ and p-. The result is shown in Fig. |l], where the mass shift of each 
contour is indicated by the number beside. As seen from the figure, the mass shift is 
not symmetric around p+ = p^ due to the second order effect. 

The opening angle could change if the decays at a distance from the origin. The 
angles between the momentum and the momentum, 6±, are given by 

= , , (5) 

l + 2(|i)cos«. + (M) 

Assume the kaon momenta are not small, so that we can approximate the kaon 
trajectories as straight lines as the kaons escape from the fireball. Using the coulomb 
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Figure 1: Shift in the invariant mass of (f) ^ K^K^ versus p+ and p_, the and 
K~ momenta in GeV/c. Numbers indicated by the contours are the corresponding 
mass shifts in MeV/c^. 



force derived from Eq. |1], 

,\ Za- f/r^ (r > R) 
\ Za-r/R^ {r<R) ' 

we obtain the momentum change in the direction perpendicular to the trajectory 
for decay outside the spherical fireball as 

^out^^r,^, b± r ^.2 , ,2\-3/2 l-COS^i 



and for decay inside as 
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where b± = r sin 6± is the distance between the straight hne trajectory and the 
origin of the fireball. Thus, the angle deflection can be obtained by 

A9± = arcsin — (9) 
P± 

and therefore the change in the opening angle by 

Ae = ef-ei = Ae+ + a^_. (lo) 

As seen from Eqs. p|-p!0|, the opening angle decreases (increases) for p_ > 
(p_ < making the mass shift more negative (positive). However, it will be 
shown in Sect. ^ that taking into account the opening angle change derived above 
does not make significant difference in the mass shift. But, the assumption of large 
kaon momenta does not always hold in the fireball rest frame. However, at low kaon 
momentum, the change in the opening angle makes larger contribution to the mass 
shift than derived from above Eqs., therefore including the opening angle change does 
not change the mass shift qualitatively. For the sake of simplicity, we will neglect it 
in the following analysis. 

3 Monte Carlo Procedure 

As mentioned above, the (f) mass shift due to coulomb effect is positive on average 
if all K^K^ pairs are accepted and reconstructed. However, the E859 spectrometer 
has a limited acceptance in which systematic decrease of the mass was reported. In 
order to implement the acceptance, we used a Monte Carlo approach. 

In the Monte Carlo analysis, we generate 0's according to exponential distribution 
in m_|_ with inverse slope 180 MeV, and Gaussian distribution in rapidity centered 
at 1.2 and with width a = 0.5 [^. The rapidity of the charge fireball in central 
Si+Au collision is yfb = 1.13, which is distinguished from the center of mass rapidity, 
Vcm = 0.97. The azimuthal angle $ of the meson is chosen randomly, and such that 
the spectrometer is around $ = 0. We assume all 0's are created at the origin of the 

m _r_ 

fireball, and decay according to e , where m, p and cr are rest mass, momentum 



5 



in the fireball rest frame and decay constant of 0. The decay angles, (O^, ^k) for 
and (tt — Ok, tt + for are selected isotropically in the (f) rest frame. 
In order to make the Lorenz transformation easy from the rest frame, {x', y', z'), 
to the fireball rest frame, {x,y,z), we select Ok = to be the direction of the 
momentum and x'-axis to be on the plane of the (p motion and the z-axis (thus 
y'-axis is parallel to the x — y plane). Given the kaon 4- momentum in the (f) rest 
frame. 
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where px is the decay momentum and rriK is the kaon rest mass, the kaon 4- 
momentum in the fireball rest frame can be obtained by 
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where f3 and 7 are the velocity and Lorentz factor of the in the fireball rest frame. 
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Figure 2: Acceptance of the E859 spectrometer at 5° (left panel) and 14° (right panel) 
with magnetic field 4A for (right contour in each plot) and K~ (left contour in 
each plot) 



We then give a ± coulomb energy shift into K"^ energy, and calculate their mo- 
menta which are measured in the spectrometer. Since the (f) was measured by E859 
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with the spectrometer at 5° and 14° and magnetic field of 4A (—4 KG) and 4B 
(+4 KG) 0, we apply only the acceptances of these settings to the resulting kaons. In 
addition, we require momenta of both kaons to be in the range of 0.4 < p < 3.5 GeV/ c 
in which the kaons are identified in the experiment 0]. We also let the kaons decay 
according to e p ^^^ where m, p and cr are rest mass, momentum in lab system 
and decay constant of kaon. In the analysis, then, we use only the 0's whose decay 
kaons are in the acceptance, within the momentum range and do not decay in the 
spectrometer before the Time-of-Flight wall which locates at 650 cm from the target 
and was used in kaon identification. Fig. ^ shows the acceptances of and K~ 
for magnetic field 4A and spectrometer angles 5° and 14°. The acceptances of 
and K~ for magnetic field 4B are identical to those of K~ and for 4A at the 
corresponding spectrometer angles. 

4 Results 

We generated 20 million (/)'s within rapidity and p± range of 0.5 < y < 2.5 and 
0.2 <p_L< 2.2 GeV/c, respectively. Out of these, about 55 K, 46 K, 10 K and 7.5 K 
0's are accepted by the spectrometer for each magnetic field polarity at 5° 2 KG, 5° 
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Figure 3: Decay angles of from {/)'s whose secondary K^K~ pairs are accepted and 
identified in the spectrometer at 14° with magnetic field 4A (solid) and 4B (dotted) 
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4 KG, 14° 2 KG and 14° 4 KG settings, respectively. As stated in Sect. ||, all 0's decay 
isotropically in its rest frame. However, decay angles selected by the spectrometer 
are strongly biased. Fig. ^ shows the decay angles of the (p^s whose secondary kaons 
are accepted in the spectrometer. 

Because of the acceptance bias, more K^K'^ pairs are populated in the negative 
(positive) (j) mass shift region for A (B) polarity setting of the magnetic field. Fig. ^ 
shows the contour plots of the accepted K^K^ pairs in the spectrometer at 14° and 
with magnetic field 4A and 4B versus and K~ momenta in the fireball rest frame. 
By comparing them to Fig. |l|, one can see that there will be negative shift in the 
invariant mass for the 4A setting, and positive shift for the 4B setting. Profiles of the 
mass shift as a function of the cj) p± are shown in Fig. ^ for 4A and 4B magnetic field 
settings with the spectrometer at 5° and 14°. 
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Figure 4: Contour plots of and K momenta of accepted pairs from (j) decays in 
the spectrometer at 14° with magnetic field 4A (left) and 4B (right) 



One might think that the positive mass shift for magnetic field B polarity should 
be larger in magnitude than the negative shift for A polarity, because one could 
just fiip the and K~ decay angles in the rest frame to have the kaons in the 
acceptance of the opposite polarity, and because the average mass shift should be 
positive. However, this is not true because fiipping the magnetic field polarity cannot 
fiip the coulomb potential of the positive fireball on or K~ . Technically, the 
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Figure 5: Average mass shift versus p^with spectrometer at 5° and 14° and magnetic 
field settings of 4A (filled-in circles) and 4B (open circles) 

argument goes as follows. The spectrometer does not accept equally a K^K^ pair 
in A setting and its image pair (with kaon momenta flipped) in B setting. Consider 
a kaon pair from a decay with momenta of (p+,p_) = (^1,^2) in the fireball rest 
frame, where pi < p2- The mass shift for this pair will be negative and the momenta 
in lab system will be {p+,P-)iab = (Pi + ^ViiVi — ^^2)- Consider another kaon 
pair from a decay with momenta of (p+,p_) = {p2,Pi) in the fireball rest frame. 
The mass shift for this pair will be positive and the momenta in lab system will be 
{p+,P-)iab = {p2 + A'p2,Pi — ^'Pi)- Because the spectrometer has the triangle-shape 
acceptance in {y, pj_) as shown in Fig. ^ it has smaller acceptance for lower momentum 
particle. Thus, the probability for the second pair to be accepted is smaller than that 
for the first pair at both polarities, therefore the average mass shift combining both 
polarities is negative. However, the effect becomes smaller when the kaon momenta 
are large, or the has large pj_. This is reflected in the mass shift at large p± at 
which the positive shift is indeed larger than the negative shift in magnitude. 

Although the average mass shift shown in Fig. |^ is large at 14° 4A setting, 
the mass shift distribution shows an asymmetric peak at a negative value close to 
zero. Fig. ^ shows the mass shift distributions with different cuts on p± of the 0. 
In addition to the 0's described earlier, the figure also includes 0's accepted in the 
spectrometer out of 10 million 0's that were generated within a more restrictive p± 
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< 0.6 GeV/c 0.6 < < 1 .2 GeV/c P^ > 1 .2 GeV/c 

Figure 6: Distributions of mass shift (MeV/c^) in three p± bins 

range of 0.2 < p± < 0.8 GeV/c. Gaussian fits to the distributions in the range of 
(—5,5) MeV/c^ indicate that the Gaussian means are significantly smaller than the 
averages. 

Including the change in the opening angle, one will get enhanced mass shift as 
discussed in Sect. 0. Fig. |^ shows the mass shift with the angle change included 
versus that without. The effect is not significant. Although the assumption of high 
kaon momenta in the fireball rest frame used in the calculation of the opening angle 




Figure 7: mass shift (MeV/c^) with opening angle change included versus that 
without 
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change does not always hold, including the opening angle change always enhance the 
mass shift, therefore does not change the conclusion qualitatively. 

Since coulomb potential is proportional to Z, the (f) mass shift increases with 
collision centrality (impact parameter). However, the increase is not significant once 
the Si+Au collision reaches a certain centrality. 

5 Limitations of the Analysis 

They are many unphysical assumptions in the analysis presented, preventing us from 
drawing quantitative conclusions. In the real world, the positive charges in the fireball 
are not static or confined in a spherical volume, and the charge density is not uniform. 
As many experimental data suggest, the fireball is expanding in both longitudinal and 
transverse directions. Therefore, the secondary kaons from decays might co-move 
with the charges longitudinally. Besides, it is not all clear whether or not the kaons 
from decays are fast enough to escape quickly from the fireball so that they are 
less effected by the transverse expansion of the fireball. To take into some of the co- 
moving effect between the charges in the fireball and the decay kaons, we alternatively 
apply an cylinder fireball with uniform charge density to the same analysis, so that 
longitudinal momenta of the kaons do not change. We assume a cylinder source with 
radius of that of the Si projectile, i?, and with length of the diameter of the Au target, 
and approximate the coulomb potential for a positive unit charge as 



where r is in cylinder coordinates. The mass shift thus obtained are more pro- 
nounced (by factor of 1.5-2) than those shown in Sect. ^ for a spherical uniform 
charge fireball. 

In the calculation presented in the previous sections, we have neglected coulomb 
interaction between the secondary kaons and the target fragments. However, we can 
estimate the effect by treating the target fragments as the source of the net positive 




(r > D/2) 
{R<r < D/2) 
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charge, thus effectively assuming all the 0's are produced at t = 0. The mass shift 
thus obtained is consistent with zero and has much less p± dependence. We note 
that this is the upper limit of the effect because any 0's produced at later times feel 
less coulomb potential from the target fragments. Although the target fragments 
have to be considered together with the fireball in a realistic calculation, they are not 
expected to make a significant contribution compared to that from the fireball alone. 

6 Prediction on Au+Au Collision 

We applied the same analysis to central Au+Au collision at 11.1 A-GeV/c. In this 
case, net charge and radius of the created fireball are Z = 158 and R = 8.2 fm, 
respectively. The fireball rapidity is yfb = 1.6, and is as same as the center of mass 
rapidity. The 0's were generated with m_|_ inverse slope 180 MeV, and Gaussian 
distribution in rapidity centered at yfb = 1.6 and with width a = 0.5. The average 
mass shift as a function of p± is shown in Fig. |] for spectrometer at 14° and magnetic 
field 4A and 4B. As expected, the shifts are about 2 times of those in central Si+Au 
collision for the same settings. 
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Figure 8: Average (p mass shift versus p^with spectrometer at 14° and magnetic field 
4A (filled-in circles) and 4B (open circles) 



In order to answer whether or not mass shift, if any, in central Au+Au collision 
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is due to non-conventional effect of chiral symmetry restoration at extreme conditions, 
one has to study differences in the data taken with various spectrometer and magnetic 
field settings. It should be pointed out that the study of the (p mass shift due to 
coulomb effect presented here does not dilute any interest in the subject of (j) mass 
shift from chiral symmetry restoration. If the chiral symmetry is restored in central 
Au+Au collision, anomalously large drop in the mass is expected and should be 
able to manifest itself from coulomb effect discussed here. 

7 Summary and Outlook 

In summary, we have reported here an analysis of mass shift due to coulomb in- 
teraction between decay kaons and a positive charge fireball. The shift is negative 
(positive) on first order of Za when the K~ momentum is larger (smaller) than the 
K^^s in the fireball rest frame. Average mass shift, summing over all possibilities 
of decay kinematics, is slightly positive due to second order effect. We implemented 
the E859 spectrometer acceptance in our calculation for settings at which systematic 
decrease in the invariant mass was observed [Q, and found that the shift caused by 
coulomb effect is in semi- quantitative agreement with the experimental observation. 
We demonstrated that the shift resulted from strongly biased decay angles selected 
by the E859 spectrometer. We applied the same analysis to central Au+Au collision, 
and predicted that the shift from coulomb effect is 2 times of that in central Si+Au 
collision. 

Due to many limitations and classic nature of the analysis, the paper offers only 
semi-quantitative look at coulomb effect on the invariant mass. To fully understand 
the effect, further analysis is clearly needed. For example, one can extend the quantum 
relativistic analysis of coulomb effect on single particle spectra in to kaon pairs 
from decays, thus analyzing the mass shift in the quantum relativistic framework. 
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